Acoustic Nature of the Boson Peak in Vitreous Silica 
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New temperature dependent inelastic x-ray (IXS) and Raman (RS) scattering data are compared 
to each other and with existing inelastic neutron scattering data in vitreous silica (v — Si02), in the 
300 -T- 1775 K region. The IXS data show collective propagating excitations up to Q = 3.5 nm~^ . 
The temperature behaviour of the excitations at Q = 1.6 nm~^ matches that of the boson peak 
found in INS and RS. This supports the acoustic origin of the excess of vibrational states giving rise 
to the boson peak in this glass. 



PACS numbers : 63.10.+a, 63.50+x, 78.70.Ck 

A peculiar characteristic of disordered solids, not 
present in their crystaUine counterpart (PhiUips, 1981), 
is an excess of low frequency states with respect to the 
predictions of the Debye theory. Inelastic neutron scat- 
tering (INS) studies on the density of states, ^(w), re- 
veal, in fact, a broad peak (the boson peak ) (BP) in the 
quantity g{i^)/uP' (Buchenau et al. 1986, Buchenau et al. 
1988, Sokolov et al. 1995) at energies 2-MO meV (15^80 
cm^^), where the Debye theory would predict a constant. 
This excess is responsible for the presence, in any glass 
studied so far, of a bump in the temperature dependence 
of Cp/T^ (Sokolov et al. 1993, Brodin et al. 1994). 

The broad boson peak feature is also observed in 
Raman scattering (RS) spectra. This can be under- 
stood considering that the first order Raman scatter- 
ing intensity, /(cj), in the harmonic approximation, is 
connected to g{uj) (Galeneer and Sen 1978): /(cj) = 
C{uj)g{uj)[n{u!, T) -f l]/w. Here n(w, T) is the Bose popu- 
lation factor and C{uj) is the photon-excitation coupling 
function, assumed to have a smooth behaviour in the con- 
sidered energy region (Fontana, Rocca and M. P. Fontana 
1987, Fontana et al. 1997). 

The nature of this excess of states is still subject to 
speculation and, at present, there are two different pre- 
vailing hypotheses (See Philosophical Magazine, special 
issue: V international workshop on disordered systems, 
Andalo, 1994). In the first one, the excess of states is ex- 
plained by the localization of the high energy vibrational 
modes induced by the static disorder in the glass. In the 
second case, collective propagating modes are thought 
to persist at high frequency, and the BP reflects their 
density of states. In the first view the BP energy corre- 
sponds to the energy of the high Q excitations, while in 
the second one the excitations are expected to propagate 
at energies above the BP energy. 

Information on the character of the excitations con- 
tributing to the boson peak can be obtained from mea- 
surements of the dynamic structure factor S{Q, oj). Using 
inelastic x-ray scattering (IXS) (Sette et al. 1995), it has 



been possible to measure S{Q,uj) of different glasses in 
the momentum-energy region of the BP, and to support 
the hypothesis of the acoustic-like character of the collec- 
tive modes and therefore of the peak itself (Masciovecchio 
et al. 1996a, Masciovecchio et al. subm.). In the specific 
case of u-silica, however, the IXS data have been alter- 
natively interpreted either with the phonon-localization 
model (Foret et al. 1997), or as propagating excitations 
(Benassi et al 1997, Masciovecchio et al. 1997). 

A better insight on the relation between the excitations 
seen in S{Q, uj) and the BP could be gained by comparing 
the temperature dependence of their energies. The basic 
idea is to utilise the observation that in vitreous silica, 
in contrast to the majority of glasses, the low frequency 
sound velocity, measured by Brillouin Light Scattering 
(BLS), (Vacher and Pelous 1976) and the boson peak 
energy (Wischnewski et al. 1998) shift towards higher 
values with increasing temperatures, and, more impor- 
tantly, the rates of change reported in the two cases are 
different. It is not the aim of this work to account for 
these peculiarities, likely to be due to the specific shape 
of the interatomic potential in silica. Rather, one can uti- 
lize the different temperature dependence of excitations 
with different Q, as seen by IXS, INS, RS and BLS, to 
establish a link between specific vibrational states and 
the features of the density of states. In particular, one 
can assess whether the excitations observed by IXS at 
energies comparable to the BP energy have either a tem- 
perature behaviour similar to that of the BP, or to that of 
the low Q excitations measured by BLS, or even another 
dependence. 

In this Letter we report temperature dependent IXS 
and RS data of vitreous silica. In the IXS spectra, prop- 
agating collective excitations are found at T w 1400 K, 
thus confirming that the high frequency longitudinal dy- 
namics has an acoustic-like nature in the whole tempera- 
ture region below melting, and up to a momentum trans- 
fer Q=3.5 nm~^, corresponding to energies more than 
twice the BP energy (E^^, w 4 meV at room tempera- 
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ture). More importantly, we find the same temperature 
shift for the BP energy and for the energy, of the 

modes at Q* = 1.6 nm-i (fi((5*) ~ E^^). This sup- 
ports the argument that the acoustic modes, probed by 
the IXS, give a significant contribution to the BP. 

The IXS experiment was carried out at the very 
high energy resolution inelastic x-ray scattering beamline 
(BL21/ID16) at the European Synchrotron Radiation Fa- 
cility. The total resolution function, measured using a 
Plexiglas scatterer at the maximum of its static structure 
factor ((3=10 nm~^), has a full width at half maximum 
[fwhrn] of 1.5 ± 0.1 meV. The Q values were selected be- 
tween 0.75 and 4 nm~^, and the Q resolution was set to 
0.2 nm~^ fwhm. Energy scans were performed by vary- 
ing the relative temperature between the monochromator 
and the analyser crystals by ±0.45 K with a step of 0.0075 
K. Each scan took about 200 min, and each Q-value was 
obtained by averaging five scans (the total integration 
time was 500 s per point). The data were normalized 
to the intensity of the incident beam. Further details on 
the IXS beamline are reported elsewhere (Masciovecchio 
et al. 1996a, b, Verbeni et al. 1996). The Raman scat- 
tering measurements were performed using a standard 
Raman laser system. Depolarized spectra were collected 
in the -300-^2000 cm~^ frequency range with 1 cm~^ fre- 
quency resolution. The 5*102 suprasil sample, purchased 
from Goodfellow, was the same 2 mm diameter rod used 
in (Benassi et al 1997, Masciovecchio et al. 1997). 

From the linear dispersion curve determined by a first 
set of IXS spectra taken at T=1375 K, it was possible 
to establish that the excitations have the same propagat- 
ing nature as that previously observed in the 300-^1000 
K region. In Fig. 1, we report as an example the IXS 
spectrum taken at Q = 1 nm~^, together with the fit- 
ting function. As in previous works, the IXS spectra 
were fitted by the experimentally determined resolution 
function convoluted with a (5-function for the elastic peak 
and a Damped Harmonic Oscillator (DHO) model (Fak 
and Dorner, 1992) for the inelastic signal. The excita- 
tion energies are reported in the inset. The data show an 
acoustic like behaviour, linearly extrapolating at small Q 
with a slope v = 6800 ± 200 m/s. This value, larger than 
the sound velocity measured by Brillouin light scattering 
(Vacher et al. 1976) implies a deformation of the acous- 
tic branch not recognized before, which will provide an 
interesting subject to further investigations. The disper- 
sion relation covers the whole BP energy region (« 5 7 
meV), and reaches energies as high as 15-;-20 meV. These 
observations, as in previous lower temperature measure- 
ments, show that, even at temperatures approaching Tg, 
it exists an acoustic-like propagating mode up to energies 
that lie well above E^p. 

A second set of IXS measurements was performed as 
a function of temperature at the fixed value Q* = 1.6 
nm^^, which was chosen since Q,(Q*) is in the range of 
Ej^p at room temperature. The IXS data and the corre- 



sponding fits are shown in Fig. 2. It is possible to observe 
directly from the row data that the intensity and energy 
of the inelastic excitation increase with T. This is em- 
phasized in the inset of Fig. 2, where the T-dependence 
of the inelastic-to-elastic intensity ratio, R{T), and the 
excitation energy, ri(Q*,T), are reported. 

To compare the T— dependence of S1(Q*,T) with that 
of other spectroscopic features, we define the scaling 
factor a,^s(r) as n{Q\T) /n{Q* ,T = OA'). Here 
f2((5*,r = QK) is the T — Q extrapolation of the mea- 
sured ri(Q*,r). The values of a-ixsC^) reported in 
Fig. 3, together with the scaling factor for the low Q 
(Q ~ 0.036 nm~^) excitation energy determined by Bril- 
louin hght scattering (Scopigno 1997), a{T)g^^. The 
difference in the T-dependence of the low- and high-Q 
excitations indicates again that the acoustic modes dis- 
persion relation is deformed from a simple linear law. 
This behaviour was hidden in the error bars in previous 
IXS and BLS studies (Benassi et al 1997, Masciovecchio 
et al. 1997). and is evidenced here thanks to increased 
sensitivity and temperature range. It is not the aim of 
the present work to discuss the origin of the observed 
deformation of the dispersion relation Q . 

This non trivial T-dependence can help us to assess 
the origin of the BP by the study of the temperature de- 
pendence of its energy. In Fig. 3 we report, in fact, the 
scaling factor a^^g (T) of the BP energy maximum, as de- 
termined from existing INS measurements (Wischnewski 
et al. 1998). The striking similarity between a^^g{T) 
and ajxg{T) provides a first strong indication that the 
states probed by IXS contribute to the BP. This similar- 
ity is further confirmed by the temperature dependence 
of the BP energy as measured by RS in the 300-;-1100 
K region. In Fig. 4 we report I{Lo)uj/C{uj)/[n{ijj,T) + 1], 
which corresponds directly to g(ijj)/uP' . We used the cou- 
pling coefficient C(aj) derived from the ratio of the INS 
(Wischnewski et al. 1998) and RS (this work) data mea- 
sured in u — Si02 at room temperature: this C(tj) is 
reported in the inset a) of Fig. 4, and it was used to scale 
all the RS spectra at different temperatures [|. The en- 
ergy of the maximum of the spectra of Fig. 4, derived 
by a local quadratic fit, is reported in inset b) and the 
scaling factors a^^{T) are reported in Fig. 3. The two 
data sets (INS and RS) are equivalent within the error 
bar. 

In presence of the deformation of the dispersion re- 



*This effect could, for example, be ascribed to the presence 
of a positive fourth order anharmonicity in the next-nearest 
neighbours interaction potentials (Scopigno 1997) or to a re- 
laxation process active in the glass. 

^It has been recently shown that the coupling function is not 
dependent on temperature in the energy region spanned by 
the BP, (A. Fontana et al. subm.). 
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lation, it is striking to observe the identical temperature 
behaviour for two spectroscopic features that in principle 
may share only a similar energy: the BP maximum and 
the acoustic excitations at Q* =1.6 nm"-'^. This demon- 
strates their common origin. Consequently, the high fre- 
quency acoustic excitations, in spite of their propagating 
nature, must have a density of states exceeding the Debye 
prediction. 

A model, that is consistent with a picture of lin- 
early dispersing propagating excitations with an excess of 
states, must take into account the effect of the topological 
disorder on the eigenvectors of the vibrational modes. It 
has been shown by molecular dynamics simulation that, 
in a glass up to the BP energy region, the eigenvector 
of a given mode j can be thought as the sum of two 
very different components, i) K plane wave like part, 
with an average momentum Qj , related to the eigenvalue 
by the linear dispersion ujj f« vQj. The Q distribution 
around Qj gets broader increasing the mode energy, ii) 
A random uncorrelated component whose spatial Fourier 
transform extends from Qj up to Q values much higher 
than the first sharp diffraction peak position with a flat 
spectral distribution (Mazzacurati, Ruocco and Sampoli 
1996). The first component accounts for the peaks in the 
dynamic structure factor and for the existence of a dis- 
persion relation. The finite width of the inelastic peaks in 
the S{Q,Lo), as that in Fig. 1, is due to the finite projec- 
tion of eigenmodes of different energy at the considered 
Q. The presence of the second component allows to sat- 
isfy the eigevectors orthogonality conditions without the 
constrain on mode counting holding for plane waves, i. e. 
without imposing the Q^ dependence of the modes den- 
sity at a given frequency. Specifically, with respect to the 
Debye behaviour, this can give rise to an accumulation 
of states at low energy. 

In conclusion, a study on the temperature dependence 
of the dynamics in vitreous silica has allowed to identify 
equivalent temperature behaviours for the maximum of 
the BP and for the high frequency continuation of the 
sound branch at the same energy. This provides argu- 
ments in favour of the acoustic-like origin for the excess 
in the density of states in vitreous silica. Reminding that 
the BP is observed basically in any disordered material, 
and that IXS has shown the presence of a propagating 
collective dynamics in all the glasses and liquids studied 
so far, it is natural to speculate on the generality of the 
conclusion reached here on the nature of the BP. 

We acknowledge A. P. Sokolov and G. Viliani for useful 
discussions. 
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FIGURE CAPTIONS 

FIG. 1 - Inelastic x-ray scattering spectrum of v — Si02 at T=1375 
K and Q = I nm~^ (o). The full line is the best fit to the 
data as discussed in the text. The dashed and dotted lines 
represent the clastic and inelastic contributions to the fit. 
The inset reports the energy of the excitations as derived 
from the fits. The dashed line is the linear extrapolation to 
small Q, and its slope corresponds to 6800 m/s. 

FIG. 2 - X-ray spectra of ti — 5102 at Q=1.6 nm~i taken at 
different temperatures. The data are shown together with 
their best fits (full line) and the individual contributions to 
the fitting function: elastic peak (dotted line), and inelastic 
components (dashed line). In the inset are reported: 1) (•) 
the inelastic-to-elastic intensity ratio, R(T)\ it follows the 
expected linear behaviour; and 2) (o) the excitation energy, 
f2(T), derived from the DHO model. The reported full lines 
are their linear fits. 

FIG. 3 - Temperature dependence of the scaling factor a(T), nor- 
malized to a(0) = 1, for the different spectra. Open symbols 
refer to incoherent measurements of the BP: INS (open cir- 
cles) and Raman scattering (open squares). Full symbols 
are the scaling factors for the excitations energies deter- 
mined from the dynamics structure factor S(Q,uj) measured 
at Q=1.6 nm-i by IXS (full circles) and at Q=0.036 nm"! 
by BLS (full squares). The dashed lines are linear guides to 
the eye to emphasise differences and similarities among the 
temperature behaviours of the four data sets. 

FIG. 4 - Examples of RS spectra taken, from to top to bottom, at 
333, 523, 823, and 1073 K. The reported data correspond to 
the reduced Raman intensities divided by the Ciu)) (Fontana 
et al. submitted) and shown in inset a). These spectra cor- 
respond to the g{oj)/uj'^. In inset b), the energy position of 
the maximum intensity at each temperature (•) is reported 
togheter with its linear fit in the low temperature region 
(dashed line). 
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Figure 1: "Acoustic Nature of the Boson Peak...", by C. Masciovecchio et al. 
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Figure 2: "Acoustic Nature of the Boson Peak...", by C. Masciovecchio et al. 
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Figure 3: "Acoustic Nature of the Boson Peak...", by C. Masciovecchio et al. 
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